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Abstract 

We examine spin vortices in ferromagnetic quantum Heisenberg models with planar 
anisotropy on two-dimensional lattices. The symmetry properties and the time evolution 
of vortices built up from spin-coherent states are studied in detail. Although these states 
show a dispersion typical for wave packets, important features of classical vortices are 
conserved. Moreover, the results on symmetry properties provide a construction scheme for 
vortex-like excitations from exact eigenstates, which have a well-controlled time evolution. 
Our approach works for arbitrary spin length both on triangular and square lattices. 
PACS numbers: 75.10.Jm, 75.10.Hk 

1 Introduction 

Vortices are a central issue in classical models for two-dimensional magnets, for a review see 
0. The dynamics of individual vortices has been studied extensively for Heisenberg models 
with easy-plane symmetry, usually combining simulations performed on discrete lattices with 
analytical approaches via continuum approximations [§, |||, ||] . These studies have led to re- 
markable insight in the dynamics of vortices in certain classical magnetic systems in terms of 
collective variables. 

With respect to statistical properties of such systems, vortices play a crucial role in the sce- 
nario of phase transitions of the Kosterlitz-Thouless type, where vortices and antivortices are 
bound in pairs below a transition temperature Tkt while they unbind above Tkt |]]- Follow- 
ing these considerations, the high-temperature phase of a planar ferromagnet was described 
by a dilute gas of topological defects, and the dynamic form factor of such a system was ob- 
tained using further reasonable approximations || . These results are in qualitative agreement 
with neutron scattering measurements on suitable quasi-two-dimensional magnetic materials. 
In particular, the dependence of the dynamic form factor on wavelength and temperature is 
found to be consistent in theory and experiment. These findings strongly support the classical 
description of such magnetic systems and in particular the existence of vortex-like excitations, 
although important aspects of this approach are still under discussion; for a critical overview 
on recent research see 0. 

On the other hand, real magnetic materials consist of quantum spins. Therefore, the question 
naturally arises whether quantum states exist which mirror the essential features of classical 
vortices. The standard answer given in the above literature and in many other papers is 
as follows: Magnetic systems with spin lengths S > 1 should be well described by classical 
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models, while for smaller spin length quantum effects become important. Nevertheless, clas- 
sical models are sometimes used also in this case, where quantum effects are approximated by 
renormalizations of coupling parameters in the Hamiltonian, see in particular |H| and refer- 
ences therein. In this work we present a concept of quantum vortices which is closely related 
to the classical limit, but takes into account the full quantum mechanics. 
The plan of this paper is as follows: In the next section we introduce the spin model we are 
dealing with and summarize some of its important properties. In section |3| we examine spin 
vortices built up from spin-coherent states. The results obtained there will lead us in section 
| to a construction of vortex-like states from eigenstates of the Hamiltonian. 



2 The model 

We consider a Heisenberg ferromagnet with planar exchange anisotropy acting on spins of 
length S on either a triangular or square lattice, 

J 



H = - - £ [SfS? + SfS] + (1 - A) Sfsf\ , (1) 
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with J > and the sum going over all pairs of nearest neighbors. We will be interested 
in planar quantum spin vortices, whose classical counterparts are known to be stable for 
sufficiently large anisotropy parameters A In the following we will always assume A = 1, 
which lies in the region of classical stabilty for both lattice types. To construct vortices we 
consider finite samples of a triangular or square lattice with open boundaries, which have 
a rotational symmetry with respect to an axis intersecting a central plaquette; examples 
are shown in figure |l|. For definiteness and brevity we concentrate on the triangular case 
throughout this paper and only briefly comment the case of the square lattice, where analogous 
results hold. 

For further reference let us briefly summarize some simple properties of such systems using 
obvious notation. The Hamiltonian is invariant under rotation of all spins around the z— 
direction in spin space, under reversal of the z-component of all spins and under appropriate 
rotations of the lattice. An adequate basis of the Hilbert space is obvious: For N spins of length 
S we define a typical eigenstate of the z-component of the total spin by \S Z ) = \Sf) 
with S z = S*. The corresponding symmetry-adapted basis vectors are given by 

\S Z , m)=M (\S Z ) + e- l 'f m R\S z ) + e~ l ^ m R 2 \S Z )) (2) 

with m £ {—1, 0, 1}, N being a normalization factor and R the operator of a clockwise rotation 
of the lattice by 2tt/3 or, equivalently, a counterclockwise cyclic permutation of the local spin 
states. The states \S z ,m) form invariant subspaces of (|l]), where the quantum numbers S z , 
m correspond to the symmetry of the model under rotations in spin space and real space, 
respectively. For S z ^ eigenstates of the Hamiltonian (|l]) with energy E are denoted by 
\S z ,m,E) and chosen to fullfill 

| - S z ,m, E) = F\S z ,m, E) , (3) 

where F = Yl i Fi is the spin flip operator which acts on each lattice site as Fi : \Sf) t— > | — Sf), 
or equivalently 

FSfF + = S^ , FS z F+ = -S z . (4) 
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Note that F is the same as a rotation of the spins by tt around the x-axis up to a possible phase 
factor; namely it holds exp(±virSf) = (±z) 2S i^. For S z = eigenstates can be characterized 
further by the spin flip symmetry and are denoted by |0, m, E, f) with / 6 {—1, 1} being the 
eigenvalue of F. Generally, each eigenstate with quantum numbers S z , m has got degenerate 
counterparts in subspaces with the same values of \S Z \, \m\. Degenerate eigenstates which 
differ only in the sign of m are related by a complex conjugation of the spin wave function. 
The case of a square lattice is obviously analogous; one simply has to infer in equation (^) 
rotations by ir/2 instead of 2ir/3 with m £ {—1, 0, 1, 2, } (mod 4). 



3 Vortices built out of spin— coherent states 

We now examine vortices which are built up from spin-coherent states on each lattice site. 
Such objects have recently been discussed by the present authors within a semiclassical ap- 
proach ||. Here we take into account the full quantum mechanics. 

In the HUbert space of a spin of length S a spin-coherent state \S;i}, tp) is defined by the 
equation 

s# >ip ■ S \S; ti,(p) = hS \S; &,<p) (5) 

for the direction s$ !(f = (sin i? cos ip, sin i9 sin ip, cost?) |L0|]. These states can be considered as 
the immediate quantum analogue to classical spin vectors. In the usual basis they can be 
expressed as 

\S;*,ip) = U{*, V ) \S) = f; ( 2S ) 2 fcos (^-)) 2S n fsin fl)Y e Mn-S) {s _ n) (6) 



n . n / V \2// V \2 

n=0 \ / 



with 



U(0, ip) = exp exp ("^) • (7) 

For our purposes we shall define here the vorticity v of a quantum state completely analogously 
to the classical case by 

v = A<P*J , Atpij = (tfj - ipi) £ [-7T, 7r[ . (8) 

The sum is taken over a closed path on the lattice in counterclockwise direction and the 
classical-like angles <pi are given by local expectation values of spin operators, 

tpi = tan" 1 (mod 2vr) . (9) 

(5f) 

Thus, the vorticity is a nonlinear functional of the underlying quantum state. In the following 
we shall restrict ourselves to the case \v\ = 1. 

We now model a planar vortex as a tensor product of spin-coherent states, 

N-l 

\i>±)= (g) \S;0i,<pi), (10) 

i=0 
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where we take = tt/2 for all i and the angles (p , to be given by the classical values as depicted 
in the examples of figure |]. This choice leads to a vortex (denoted by with v = 1), while 
the mapping </?j — ^ converts it into an antivortex \ip-) with v = —1. From (Q) one sees 
that and |^_) are related via a complex conjugation of the spin wave function, which is 
the same here as a spin flip 

\i, + ) = F\^). (11) 



Before examining further the quantum states ( 10 ) let us briefly remark on the classical vortex. 



As seen in figure |], in the small system a) all directions of the classical spins can be derived 
by intuitive symmetry arguments and are the same as in the well-known continuum solution 
<p(x, y) = tan _1 (y/x) + constant with x, y denoting cartesian coordinates in the plane. In the 
system b) the same holds for the inner lattice sites labelled by 0, 1, 2 and 5, 8, 11, but not 
for the outer ones. E. g. it is easy to see that the sum of the classical vectors on 3 and 4 
must be always parallel to the spin on 0, but the exact value of, say, (j)^ must be calculated 
in detail and turns out to be different from the continuum solution. In system c) of figure [l] 
one has the same situation for the sites 4, 5 and so on. Note also that the classical vortex is 
a static solution only if its center coincides with the center of the system, because otherwise 
its image vortices created by the boundaries cause a movement of the vortex center. Within 
the continuum approximation of the system this situation is completely analogous to two- 
dimensional electrostatics. 

Now we analyze the states \ip±) with respect to the symmetries of the Hamiltonian. Let us 
concentrate again on the triangular case. The scalar product of the vortex (|l0|) with a typical 
basis vector (0) is 

(S z , rr#±) = M ((S z \4> ± ) + e + ^ m (S z \R + \4>±) + e + ^ m (S z \(R + ) 2 ^ ± )) . (12) 

The application of R + = R^ 1 on \ip±) is the same as a counterclockwise (clockwise) rotation 
of each local spin-coherent state by an angle of 27r/3, i. e. all angles <pi in ( |i~0| ) get a turn of 
±27r/3. Therefore, with the help of (||) one finds for integer spin length S 

(cx{S*\4>±) S z Tm = (mod 3) 
(5 ' m|Vj±) { =0 otherwise (13) 

and similarly for half-integer S 

(S «,„ #± }( K<S1fa> s 2 + f=F-» = ° (14) 

x ' ir±/ [ =0 otherwise v ; 

These relations determine the invariant subspaces of the Hamiltonian in which the vortex state 
|?/>±) has non-vanishing overlap. We therefore call them selection rules. To cover the case of 
the square lattice one simply has to replace (mod 3) with (mod 4). 

The square moduli of the coefficients in the expansion (||) form a binomial distribution of 
range 25 with a probability parameter p = sin 2 (?9/2). This is the probability distribution for 
the results of measurements of the z-component of an individual spin being in the state (^). 
By an elementary theorem of stochastics, the distribution of a finite sum of quantities which 
are binomial-distributed with a common parameter p is again of the binomial type with the 
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same parameter and a range just given by the sum of the individual ranges. In a planar vortex 
we have pi = sin 2 (-7r/4) = 1/2 for all i and therefore 



££|^m,^)|>=^j^. (15) 

The sum goes over all eigenstates having S z as the quantum number of the total spin; for 
S z = the states |0, m, E, f) have to be inserted and summed over / as well. The mean value 
of this symmetric distribution is of course zero and the square variance is given by 

CAT 

(AS Z ) 2 = Y . (16) 

According to the central limit theorem of stochastics, the expression ( |T5| ) approaches a Gaus- 
sian shape for large SN, where 

lim £ £|<S z ,m,£#±>| 2 = fdxg{x) (17) 

°° aAS*<S*<bAS* m,E Ja 

with the normalized Gaussian distribution g{x) = l/\/2n exp(— x 2 /2) and real numbers a < b. 
If one fixes a certain value of m in the above summations, only every third value of S z gives 



a non-vanishing contribution because of the selection rules (13), (|14|). Thus we find 

/+oo 1 
dxg(3x) = ~, (18) 

i. e. for an infinite system, N — > oo, or in the classical limit, 

7l — ^ , S — > oo , ftS 1 = constant , (19) 

the (anti-) vortex \ip±) has the same square overlap in all subspaces characterized by different 
rotational quantum numbers m € { — 1,0, 1}. 

The same arguments hold for the square lattice with 3 to be replaced with 4 in fllq). 
In summary, the above equations characterize the states \ip±) with respect to the symmetries 
of the system. In figure § we have illustrated the results for the system shown in figure [j]a) 
(N = 6) and S = 5/2. 

Next we analyze the states \ip±) with respect to the exact eigenstates of the model (Q). To 
this end we have numerically diagonalized the full Hamiltonian for small systems, i. e. have 
computed all eigenvalues and eigenvectors in the invariant subspaces. This procedure can be 
done with today's computers for the system a) in figure [l] for spin lengths S = 1/2, 1,... ,5/2, 
while for larger lattices like b) and c) one is still restricted to S = 1/2. 

Let us first consider the system shown in figure |l]a) . The expectation value of the Hamiltonian 
is 

W±m±) = ~j(nsf. (20) 

Its variance has been obtained in reference |§ and reads here 

AH = J(hS) 2 ^. (21) 
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As it must be, this quantity vanishes in the classical limit (|T9|). Figure |3| shows histograms of 
the square overlap of \ip±) with the eigenstates of the Hamiltonian and the density of states 
for S = 2 as a function of the energy; for other spin length S = 1/2, 1, . . . , 5/2 the data looks 
qualitatively similar. The time evolution of the system being initially in the state \ip±) can be 
followed in terms of 

<V±(0)|<Mt)> = (^±|e-i H *|V>±), (22) 

which is essentially the Fourier transform of the data shown in the upper diagram of figure ||[ 
Therefore this quantity decays on a time scale given by the uncertainty relation 

AH At > ^ , (23) 

which is in full agreement with our numerical findings even for comparatively large spin lengths 
like S = 5/2, where a classical-like behaviour of spin systems is often assumed. Thus, even for 
large S the time dependence of the scalar product ( |22] ) is the same as for any usual dispersive 
wave packet, in contrast to the classical vortex which is a nonlinear coherent excitation. This 
is due to the fact that the classical limit (^) is not approached properly by taking a large 
spin length but keeping fi as finite as it is. 

Inspite of this general statement, the analysis of local expectation values 

(Sf (t)> := ^M Ht Sfe-^ m \^ ± ) (24) 

shows that certain features of the initial vortex structure nevertheless remain present in the 
time evolution of the state. First note that the local spin expectation values on sites i, j 
connected by a rotation of the lattice, i. e. = RSfR + , are related by 

(Sffl) = (Sf(t)) , (S+(t)) = e**(S?(t)) , (25) 



as it is intuitively obvious and can be shown easily with the help of (|l3|), (14). Therefore, the 
vorticity of the central plaquette is necessarily conserved. 

Let us now discuss the time evolution of the vortex in more detail. For the system in figure 
|l]a) our numerical results are as follows. The expectation values (Sf(t)} are strictly zero for 
all times, lattice sites and spin lengths S = 1/2,1, ... ,5/2. This is surprising since only the 
^-component of the total spin is conserved due to symmetry. In figure || we have plotted 
the time evolution of the in-plane spin components for the state and S = 2. The upper 
diagram shows which may be seen as an 'effective spin length'. This quantity decays 

for both classes of lattice sites on a time scale given by (^) to comparatively small numbers 
and even becomes zero for certain times. The lower diagram shows the direction angles ipi 
calculated from @. Surprisingly, these angles remain constant up to changes of ir, which 
occur when |(>S^"(i)}| goes through zero, i. e. the spin expressed by its expectation values 
reverses its direction. The times when such reversals occur are not identical for both classes 
of sites, but apparently strongly correlated. Note that this gives rise to quantum fluctuations 
of the vorticity as defined in (p|). For instance, if the inner lattice sites 0, 1, 2 have undergone 
such a reversal while the outer ones have not, the vorticity on the three outer plaquettes is 
changed from to —1, while the vorticity of the central plaquette is preserved as mentioned 
before. An evaluation of the first 1000 time units after starting the dynamics shows that 
in about 80% of this time intervall the vorticities on all plaquettes have their initial values, 
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while in the remaining time the vorticities of the outer plaquettes are changed to ±1. This 
shows the strong correlation in the spin dynamics also seen in figure ||. It is an interesting 
speculation whether such fluctuation phenomena are related to the sponteneous appearance 
of vortex-antivortex pairs (in larger systems), which is well-known from classical spin models. 
The findings described above hold similarly for all spin lengths S = 1/2, 1, . . . , 5/2 and both 
types of states \ip±)- 

In the system of figure |l|b) some new observations are made. As already mentioned the 
numerical analysis is restricted to S = 1/2. The spins on the inner lattice sites 0, 1, 2 and 5, 
8, 11 show completely the same behaviour as in the system described before, while the time 
evolution of spins on the outer sites, say 3 and 4, is different. Here small z-components (Sf(t)) 
arise which are plotted in the upper diagram of figure ||. For symmetry reasons these quantities 
differ in sign on sites which are inequivalent under rotation, since the expectation value of the 
z-component ot the total spin is constantly zero. Moreover, also the in-plane angles ipi are 
not conserved (up to changes by tt) as shown in the lower diagram. But remarkably, the 
sum (Ss(t)) + (iS4(t)) is always parallel or antiparallel to (So(t)) with the orientations being 
correlated in a similar way as described before. This is also a strong reminiscence of the 
classical vortex structure. 

The differences in the behaviour of the spins found in the system of figure |l|b) are a surprising 
parallel to our previous remark on spin directions in the classical vortex. Here the outer lattice 
sites are also distinguished from the inner ones, since their spins are not described by the static 
continuum solution. 

Summarizing, we have demonstrated that the states \ip±), although they show dispersion, 
preserve important properties of classical vortices. 



4 Vortices constructed from exact eigenstates 

Since the time dependence of the spin vortices presented in the last section is rather com- 
plicated, it is desirable to find vortex-like quantum states which have a well-controlled time 
evolution. To this end the symmetry rules fll3|), (14) provide a simple construction scheme. 
It is useful to distinguish three different cases: 

(i) Triangular lattice, SN integer: A vortex-like quantum state is given in terms of exact 
eigenstates of the Hamiltonian by the following ansatz: 



f (a_i| -l,Tl,#i) +a |0,0,£o,/> + ai|l,±l,£i)) 



|x±) = ,> i ,, i i, <*-i| - 1, Tl, E x ) + a |0, 0, E J) + ai |l, ±1, Ex) (26) 

^/|a_i| / + |aor + |ati| 

This is a linear combination of eigenstates which has nonzero amplitudes only for quantum 
numbers 'allowed' by the rules (|i"3|) , ( |l4|) and is restricted to the most important contribu- 
tions with 15*1 < 1. From each invariant subspace only one eigenstate is involved. The 
states with S z = ±1 are chosen to be degenerate (which is always possible). Thus ( p6[) 
is effectively a two-level-system with an internal frequency u = {E\ — Eq)/H. Denoting 
\x±(t)) = expC - (i/fi')Ht)\x±) one finds similarly as in (|25f) : 

(x±(t)\Sl\ X ±(t)) =e ± ^( X ±(t)\St\x±(t)), (27) 
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Table 1: Expectation values of in-plane spin components in the state \x+) constructed from 
elementary excitations (see text) 

where the lattice sites i, j are related by a rotation. Therefore the central plaquette carries a 
constant vorticity of v = ±1. For = |ai| it holds: 

<X±(*)l^lx±(t)> = (28) 

for all sites i and the expectation values of the in-plane components are given by 

2|a 1 ao|(0,0,£ ,l|S+|l,±l,£ 1 ) 
2|ai| 2 + |«o | 2 

• exp Q (0_i - +1 )J cos (ut + (j) - ~ (0_i + + i)J , (29) 

where we have inferred = |q^| exp(z^) and assumed / = 1 for simplicity; the case / = — 1 
leads only to unimportant additional phase factors. To derive (p8|), (pS|) the relations (|3|), (||) 
have been used. 

Thus, the vector of the local expectation values of the spin components has a constant di- 
rection (up to reversal) on each site while its length varies harmonically with the frequency 
u in time. Differently from the vortices constructed from spin-coherent states, all vectors of 
expectation values lie strictly in the plane and their reversals, i. e. the zeros of their length, 
occur simultaneously on all lattice sites. Therefore no fluctuations of vorticity arise. 
The above construction works on triangular lattices of the type shown in figure |l| and of arbi- 
trary size. It provides quantum states with a very simple time evolution and typical properties 
of vortices. 

To illustrate this, let us return to the system la). Figure [6| shows the lower part of the spec- 
trum for S = 2 as a function of S z . The ground state has quantum numbers S z = 0, / = 
and is part of a band of states with m = 0. Well separated from this we have a band of 
degenerate states with m = ±1 and next a more or less continuum-like distribution of states. 
This turns out to be qualitatively the same for the other spin lengths considered here. The 
choice of states to be used in ( p6| ) is certainly not unique. To give a definite example, let us 
choose the combination with lowest possible energy, i. e. we take the state with S z = 0, m = 
to be the ground state and the other ones from the excited elementary band. Note that the 
expectation value of energy for such a linear combination is much lower than the expression 
( pO] ) , which corresponds directly to the energy of a classical vortex. In table |l] we present the 
expectation values of the in-plane spin components for different spin lengths in the state \x+) 
at time t = 0, where we have set all \ai\ = 1 and adjusted the phases 4>i in a manner that 
the argument of the cosine in (^) vanishes and the spin on the site has (fo = 0. Obviously, 



(x±(t)\st\x±(t)) = 



8 



the 'effective spin lengths' are strongly reduced compared with the original ones. This was 
also found in the previous section in the time evolution of a vortex built out of spin-coherent 
states (cf. figure The importance of this effect has also been stressed by other authors 
using a variational approach to the spin dynamics [11, 12 1. 

As mentioned above, the central plaquette has vorticity v = 1, and the directions on the other 
sites also exactly mirror the classical vortex structure. This is a non-trivial property since 
the only strict relation between these directions is given by (|27|). This observation can also be 
made for other choices of eigenstates, mainly from the lower part of the spectrum, but for an 
arbitrary linear combination of the form (|26| ) this is not the case. 

Thus we have demonstrated the existence of vortex-like quantum states built up from ele- 
mentary excitations , whose energy is much lower than the semiclassical vortex examined in 
the previous section. The energy of a single (semi-) classical vortex is known to grow logarith- 
mically with the size of the system. The construction presented here relies only on symmetry 
properties and works for arbitrary system size. Therefore, the energy of the vortex-like quan- 
tum state discussed in the above example must be assumed to remain in the order of the 
exchange integral J (or at least finite) even in an infinite system, since otherwise the energy 
difference between the lowest and the first excited band of quantum states would have to grow 
with increasing system size to macroscopic values (or infinity), which is completely unlikely. 



We only sketch the remaining two cases . 



(ii) Triangular lattice, SN half-integer: Here all values of S z are also half-integer and a 
vortex-like quantum state can be constructed as (cf. (|i~4])): 

\X±)= A I , | j^ fa-il-i ^,E) + a 1 \\, ±1,E)) , (30) 
V|a-i| + \ a i\ \ 2 1 / 

which also has the properties (p7|), ( p8|) for |a_i| = |«i|. Therefore the spin structure ex- 
pressed in local expectation values is also planar and the central plaquette carries a vorticity 
of ±1. The spin structure on other plaquettes depends on further details and can be examined 
as above. As the involved eigenstates with different quantum numbers are chosen degenerate, 
we obtain an exact eigenstate which has typical features of a vortex, at least with respect to 
its center. 



(Hi) Square lattice, SN necessarily integer: This case is merely analogous to (i) with the 
extension that the eigenstate in ( p6| ) with S z = may may be chosen from two different 
subspaces (m = or m = 2). This slightly generalizes the selection rules (|l3|), Jl4|), which 
give only one of these two possibilities, and leads also to different vorticities. 



5 Conclusions 

In this work we have examined planar quantum spin vortices in ferromagnetic Heisenberg 
models taking into account the full quantum mechanics. 

Vortices built up from spin-coherent states are studied in detail. These objects can be seen 
as the immediate quantum analogue of a classical static vortex on a discrete lattice. Impor- 
tant results on their symmetry properties are given by the relations (13), (0), (|ilt ) and are 
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illustrated in figure ||. The time evolution of such vortices is in general quite complicated 
and, from a global point of view, typical for quantum mechanical wave packets. On the other 
hand, a detailed numerical study of the local spin expectation values shows that important 
properties of the initial classical-like vortex structure are conserved. This may be viewed 
as a reminiscence of the topological character of the classical vortex solution, although such 
topological arguments do not apply strictly in a discrete system. 

Our symmetry analysis heavily relies on the symmetry of the underlying lattice sample with 
respect to the vortex center as shown in figure |l[ To characterize a vortex whose center lies 
not in the center of system, one may consider a subsystem which has this property. The 
quantum vortex state projected onto the Hilbert space of this subsystem should have similar 
properties as found here, e. g. the selection rules (|l3|), ( |l4| ) should hold approximately, but 
not exactly. We expect the deviations from these selection rules to result in a movement of 
the vortex center, as it is well-known from the classical vortex. 

Our observations should generally raise the confidence in applying the classical Kosterlitz- 
Thouless theory to real magnetic systems (consisting of quantum spins) in the spirit of an 
effective field theory and with respect to critical behavior, where many details of the system 
can be expected to be unimportant. 

The symmetry properties of such vortices lead to a 'reduced' construction of vortex-like exci- 
tations in terms of exact eigenstates of the Hamiltonian as described in the foregoing section. 
We obtain vortex-like quantum states involving only two different energy levels or, in partic- 
ular cases, exact eigenstates having vortex-like features. Moreover, we find such states, which 
have unambigously the properties of a classical vortex and a very simple time evolution, even 
at energies which are much lower than the classical vortex energy. This may, on the other 
hand, indicate some important modification of the role of vortices in the quantum system 
compared with the classical case. Concerning this issue further investigations are desirable, 
which may possibly start from the construction scheme of vortex-like excitations in terms of 
exact eigenstates given in this work. 

The approach presented here is expected to be useful also for other cases like non-planar 
vortices in ferromagnets or vortices in antiferromagnets. 
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Figure 1: Classical vortex structrures on lattice samples which have a rotational symmetry 
axis intersecting the center of the vortex 
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Figure 2: Total square overlap of \ip + ) in invariant subspaces of the Hamiltonian. The dashed 
line is a Gaussian with width AS" 2 . 
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Figure 3: Square overlap of \ip±) with eigenstates of the Hamiltonian (upper diagram) and 
density of states as a function of the energy for the system of figure la) and S = 2 
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Figure 4: Time evolution of local expectation values in the state for the system of figure 
la) and S = 2 
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Figure 5: Time evolution of local expectation values in the state for the system of figure 
lb) and S = 1/2 
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Figure 6: The low-lying part of the spectrum of the system of figure la) for S = 2 
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